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Abstract 

Thin liquid sheet flows are excellent candidates for 
use in space radiator systems. Surface tension forces at 
the edges of a thin liquid (-200 p.m) sheet flow result in a 
triangularly shaped sheet. Such a geometry is ideal for an 
external flow radiator. Since the fluid must have very 
low vapor pressure, Dow Coming 705 silicone oil was 
used and the emittance of a flowing sheet of oil was 
determined by two methods. The emittance was derived 
as a function of the temperature drop between the top of 
the sheet and the coalescence point of the sheet, the sink 
temperature, the volumetric flow and the length of the 
sheet. The emittance for the oil was also calculated using 
an extinction coefficient determined from spectral trans- 
mittance data of the oil. The oil’s emittance ranges were 
high, ranging from 0.67 to 0.87 depending on the sheet 
thickness and sheet temperature. The emittance derived 
from the temperature drop was slightly less than the 
emittance calculated from transmittance data. 

Introduction 

The liquid sheet radiator (LSR) is an external flow 
radiator similar in several ways to other external flow 
radiators such as the liquid droplet (LDR) 1 and liquid 
belt radiators. 2 All of these radiator concepts potentially 
have lower mass than solid wall radiators such as pumped 
loop and heat pipe radiators. They are also nearly im- 
mune to micrometeoroid penetration. However, the LSR 
has the added advantage of simplicity. Surface tension 
causes a thin (100 to 300 p,m) liquid sheet to coalesce to 
a point (Fig. 1). As a result, the sheet flow has a triangu- 
lar shape. Such a triangular sheet is desirable since it 
allows collection of the flow at a single point. To obtain 
a similar triangular area for the LDR requires very accu- 


rate aiming of many (10 5 to 10 6 ) droplet streams. Sim- 
plicity of the LSR should also result in lower mass for 
the LSR compared to the LDR and liquid belt radiators. 

A major problem for all external flow radiators is 
the requirement that the working fluid be of veiy low 
(~10 8 torr) vapor pressure in order to keep evaporative 
losses low. As a result, working fluids are limited to cer- 
tain oils (such as used in diffusion pumps) for low tem- 
peratures (300 to 400 K) and liquid metals for higher 
temperatures. 

Previous research on the LSR has been directed at 
understanding the fluid mechanics of thin sheet flows 3,5 
and assessing the stability of such flows, especially with 
regard to the formation of holes in the sheet. 6 The latest 
topic investigated was the emittance of thin sheet flows. 
The emittance was calculated from spectral transmit- 
tance data for the Dow Coming 705 silicone oil. The 
emittance was also derived as a function of measurable 
quantities, namely, the temperature drop between the top 
of the sheet and the temperature at the coalescence point 
of the sheet, the sink temperature, the volumetric flow 
rate and the length of the sheet. Since the temperature 
drop was on the order of a few tenths of a degree, very 
accurate measurement was needed. In addition, because 
the temperature drop was so small, temperature fluctua- 
tions within the oil were a potential problem and this was 
investigated. 

Analysis 

Theoretical Sheet Emittance from Transmittance Data 

The following expression is given for the spectral 
emittance, of an infinite sheet, 3 

e X = 1 - T X =1 - 2E 3K T ) W 
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where %x is the spectral transmittance, is the extinc- 
tion coefficient, x is the thickness of the sheet and E 3 (x) 
is the exponential integral. 

E,(x)= f u e (_x/u) du (2) 

J J 0 

The total hemispherical emittance is defined as follows, 4 

j^e, <X.T)d>. 

c — *0 ^ 


surface and it is assumed to behave as a black body 

(e~=l). 

Before proceeding, a number of simplifying assump- 
tions can be made. The Froude number (the ratio of the 
kinetic energy to gravitational potential energy) for this 
experiment is high enough that gravitational effects on 
the sheet geometry can be neglected. In that case, the 
sheet shape is exactly triangular. 5 If heat transfer due to 
conduction or convection is neglected (this is reasonable 
provided a sufficient vacuum is maintained), then in the 
steady state, the radiated power must balance the en- 
thalpy change, 


where e^b is the black body hemispherical spectral emis- 
sive power, 4 T is the temperature, X is the wavelength, 
and a s b is the Stefan-Boltzmann constant (5.67XKT 8 W/ 
m 2 K 4 ). 

The black body hemispherical emissive power is 
expressed as, 


e 
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Here C, = hc 0 2 = 0.59544xl0 8 W |im 4 /m 2 and C 2 = he,,/ 
k = 14388 pm K; where h is Planck’s constant, c 0 is the 
speed of light and k is the Boltzmann constant. 

From the measured spectral transmission, Tx, the 
extinction coefficient, a^, is calculated using Eq. (1). 
The hemisperical emittance, £^, is then calculated for 
any thickness, x, and any temperature, T, using Eqs. (1) 
to (3). 


Prad=PQC p ( T l- T 2) (6) 

Here p is the density of the fluid, C p is the specific heat, 
and Q is the volumetric flow rate. If it is assumed that the 
emittance is a constant for the sheet and that T varies 
only in the z-direction then Eqs. (5) and (6) yield the 
following: 

aT=T '- T > = ^JoV- T -)H) dz <7) 

Integration of Eq. (7) requires knowing T(z), which 
is found by solving the energy equation. Assuming that 
temperature gradients exist only in the flow direction 
yields the following for the steady state energy equation: 

PQCpdT = -2a sb e(T 4 - T 4 ) (8) 

Wdz sb V “/ 


Experimental Sheet Emittance 

The flow of the liquid sheet radiator is depicted 
graphically in Fig. 1 . The working fluid flows with ve- 
locity, w G , from a slit of prescribed width, W, and thick- 
ness, X, at one end of the device (z = 0) to form a 
triangular sheet of liquid. 

The planform area of the sheet acts as a radiating 
surface. When surrounded by a sink temperature, T^, 
lower than the initial fluid temperature, the fluid tem- 
perature decreases as a function of the vertical position, 
z, until reaching the coalescence point. The power radi- 
ated from the sheet is given by the following. 


By defining the following dimensionless parameters: 


h PsbEWLT, 3 

PQCp 


and 
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Eqs. (7) and (8) can now be rewritten as follows, 


r n d=2'’,bJ A E ( T4 - T -) dA <5) 

where T is the temperature of the sheet, A is the sheet 
area and £ is the total hemispherical emittance. The 
factor of two arises because both sides of the sheet radi- 
ate. The background in the experiment is a black painted 
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For 0 < < Ti, the solution to Eq. (12) subject to the 

initial condition T = Ti at z = 0 is given by: 
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(T/Tj - T^/Tj )(l + T^/Tj ) 
(T/Tj + T^/Tj )(l - T^/Tj ) 
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Equation (13), which gives T/Ti as a function z, is 
used to carry out the integration in Eq. (11). 

The solution for the temperature drop, AT, can be 
solved in the limit of low background temperature 
(T^ -» 0). The method outlined in Appendix B was used 
to solve Eq. (1 1) for AT in the limit of low background 
temperature (T** -> 0). This result is given by Eq. (B-4): 


AT 4h + l-(6h + l) 2/3 

lim = ^ — 

T^oh^ 4h 2 
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compares the exact solution (obtained by a numerical 
solution of Eqs. (11) and (13) using Newton’s method) 
with the constant temperature approximation, Eq. (14), 
at a worst case value of h = 0.004, corresponding ap- 
proximately to Ti = 400 K, W = 23.5 cm and x = 100 pm. 
This represents the highest value of h that is encountered 
in the experiments. As Fig. 2 shows there is excellent 
agreement between the exact and approximate solutions. 

The percent deviation between the exact and 
approximate solutions has an upper bound given by 
Eq. (B-14): 


( AT )approx ( AT )exact^8 h 
( AT >approx 3 

For the condition h = 0.004 considered in Fig. 2, the 
deviation is approximately 1 percent. Thus, the constant- 
temperature approximation is sufficiently accurate when 
calculating the temperature drop. 

If a measurement for AT is available then the emit- 
tance can be solved for by a straight-forward integration 
of Eq. (8) for constant T = Tp 


The range of h appropriate for test conditions is 
0.0004 < h < 0.004. For these veiy small values of h, 
Eq. (B-4) is found to be, 
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Thus, AT « hTi which is the same solution obtained 
from Eq. (1 1) if T « Ti and T** — » 0. Therefore, in the 
case To* — > 0, it is a good approximation to assume T ~ 
Tj in the integration of Eq. (1 1). The relative error in this 
constant temperature approximation was found numeri- 
cally in Appendix B to be proportional to h, 

exact 

tJSo (AT) approx 

= 4h 2 — 4h- 1 4(6h + 1 f 3 = 8 h 
4h 2 3 

In experimental conditions, T** is not likely to approach 
zero, it may be closer to Tp An analysis of T** T\ is 
given in Appendix B. It was found that the constant 
temperature approximation is still valid. The limit on the 
error in making the approximation was <8h/3. Since h 
encountered in the experiments is 0.004 to 0.0004, the 
constant temperature approximation is valid. Figure 2 


The emittance is calculated directly from eq. (14) 
provided measurements are available for the three tem- 
peratures Ti, T 2 , and To*; the flow rate, Q; and the sheet 
length, L. The specific heat, C p , and the density, p, are 
known as functions of temperature. 

Results 


Theoretical Sheet Emittance from Transmittance Data 

Transmittance data was taken for a sample of the 
Dow Coming 705 silicone oil used in the emittance 
measurement. The measuring device was an FTIR spec- 
trophotometer. The spectrophotometer measures the over- 
all transmission of the oil plus the transmission of the 
windows containing the oil (Fig. 3). Therefore the potas- 
sium bromide (KBr) window transmittance and reflec- 
tance must be taken into account in determining the 
extinction coefficient. 

The oil transmittance is shown as a function of 
wavelength in Fig. 4 and the window’s transmittance is 
shown as a functions of wavenumbers in Fig. 5. From 
this data, the extinction coefficient was calculated ac- 
cording to Eq. (A-25). It is shown in appendix A that the 
extinction coefficient needed for the spectral emittance 
can be expressed as, 
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Tx is the overall spectral transmission through the 
oil and the two windows, r aw is the reflectance at the air- 
window interface, r mw is the reflectance at the media- 
window interface (in this case the media is silicone oil), 
T w is the internal spectral transmission of the windows 
and d is the thickness of the media. 

For the temperature range (300 to 400 K) where the 
silicone oils can be used, exb is negligible for X < 2.5 
and also for X > 70 pm. Because of this e was calculated 
from ax for the wavelength region 2.5 < X < 70 pm. The 
extinction coefficient is plotted as a function of wave- 
length in Fig. 6. The hemispherical emittance was then 
calculated using Eqs. (1) to (4) for various temperatures 
and sheet thicknesses. The extinction coefficient has a 
maximum around 9 to 10 pm. For 300 to 400 K, the 
black body hemispherical emissive power is a maximum 
at 8 to 10 pm. Since the extinction coefficient is a maxi- 
mum in this wavelength region, the total hemispherical 
emittance will be large for the 300 to 400 K temperature 
region. 

Experimental Sheet Emittance 

A schematic of the test facility that was used is 
shown in Fig. 5. It consists of a 30 cm inner diameter 
stainless steel pipe 3.5 m long. The axis of the pipe is 
aligned with the gravity field. Vacuum conditions exist 
within the pipe with the pressure about 0.02 to 0.04 torr. 
At these conditions, aerodynamic drag on the sheet flow, 
as well as heat transfer due to conduction and convec- 
tion, is negligible. Dow Coming 705 silicone oil is 
pumped up to a plenum above the slit. Within the ple- 
num, vibrations and temperature fluctuations are damped 
out. The design of the slit plate is shown in Fig. 6. The 
fluid temperature was maintained at 373 K, while the 
sink temperature was maintained with a coolant at 
293 K. Four different slit plate sizes were tested: 


100 pm x 12.5 cm 
150 pm x 23.5 cm 
200 pm x 23.5 cm 
300 pmx 18.75 cm 

With these conditions, the temperature drop is very 
small, around 0.3 K. For this reason, carefully calibrated 
ultrastable thermistors were used to measure T i and T 2 . 
The thermistors are responsive to 0.01 K temperature 
changes. The upper thermistor was located directly above 
the slit in the plenum. The lower thermistor was located 
on an actuator probe which swung around horizontally 
and travelled vertically to position the thermistor at the 
coalescence point. The sink temperature does not need 
the accuracy that the temperature drop needs; thus, to 
measure the sink temperature, thermocouples were posi- 
tioned on the wall of the vessel and their readings were 
averaged to a single sink temperature. To measure the 
volumetric flow rate, Q, the pressure drop across a cali- 
brated orifice in the oil supply line was measured. Q is 
calculated according to the equation 


Q = V^P (19) 

The orifice has been calibrated with Dow Coming 
705 silicone oil at a temperature of 358 K for which 
y = 2.92 cm 6 /s 2 Pa. The density and specific heat of the 
oil is also needed to calculate the emittance. Best fit 
equations were made using data from Dow Coming 
(Ref. 9). These equations are, 

p = [l569 — (1.625T) -1 j kg/m 3 (20) 


C p = [-9017.3 + 69.888T - .15274T 2 

+ 1.1157 X 10~ 4 T 3 jj/kg K (21) 

There are two ways to determine the length of the 
sheet. In Ref. 5, it was shown that in the case where 
gravity can be neglected, 



where We is the Weber number, the ratio of the dynamic 
pressure to the surface tension pressure. 


We = 


P^O 

o 


(23) 
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where w Q is the velocity of the oil along the length of the 
sheet. It can be found by dividing the volumetric flow 
rate by the initial width and thickness. The initial thick- 
ness is x 0 and c is the surface tension. A best fit equation 
of the surface tension was made from Dow Coming data 
and is expressed as, 

a = [.08614 - 1 .5379 x 10 T]N/m (24) 

The length can also be measured directly. The actuator 
probe that the lower thermistor is located on records the 
vertical position of the coalescence point. 

Measuring the small temperature drop proved to be 
the most difficult part of the experiment. Although the 
lower thermistor was located on an actuator probe and 
had some range of motion, it was difficult to position it 
in the flow. Furthermore, at the coalescence point where 
the thermistor was located, the flow was often quite 
erratic. Fortunately, it was evident when the thermistor 
was no longer immersed in the flow because the tem- 
perature drop would increase considerably. This large 
increase was due to the lower thermistor radiating to the 
cold wall of the containing vessel instead of measuring 
the temperature of the oil. If the experiment were to be 
continued, an improved design for the measurement of 
the lower temperature would be implemented such as a 
cup with a hole in the bottom that collects the oil, with 
the thermistor located beneath the hole. Thus, the ther- 
mistor would always be in the flow. It was also a concern 
that temperature fluctuations upstream of the slit plate 
might cause error. 

An experimental investigation to measure the typical 
fluctuations was also carried out. The thermistors re- 
sponse time is 1.5 seconds, thus, only the temperature 
fluctuations with a period greater than 1 . 5 seconds could 
be analyzed. Experiments were conducted where Ti 
was measured at 1.5 second intervals. The average tem- 
perature change in this interval was found to be 0.031 K 
and the largest temperature change was found to be 
0.054 K. Fluctuations of 0.05 K or greater occurred 
17 percent of the time. 0.05 K is large for our purposes, 
however, if the period of the fluctuations is much greater 
than the flow time (the time it takes a particle to tranverse 
the sheet or L/w 0 ), then the temperature fluctuation’s 
effect on AT will not be significant. The flow time was 
measured and found to be 1.54 msec. Since the fluctua- 
tions within the system are most likely generated from 
the heaters cycling on and off, the fluctuations most 
likely have periods much greater than 1.54 msec. Thus, 
it was assumed that temperature fluctuations will not 
affect the temperature drop due to radiation. 


The resulting experimental emittance using Eq. (13) 
along with theoretical emittance is shown in Fig. 9. The 
emittance of the oil in the 300 to 400 K region proved to 
be quite high. The experimental emittance of the sheet is 
between 0.74 and 0.85, depending on the sheet thickness 
and oil temperature. The emittance calculated from the 
transmittance data agreed closely with the experimental 
values. At T = 373 K there was a difference of 0.02 in the 
emittance measured for the 100 pm sheet, a difference of 
0.01 in the emittance measured for the 150 pm and 
200 pm sheet, and a difference of 0.015 in the emittance 
found for the 300 pm sheet. The experimental values 
were slightly less than the theoretical values for all sheet 
thicknesses. As shown in Fig. 9, the emittance increases 
with an increase in temperature and an increase in sheet 
thickness. As the temperature increases the maximum 
black body hemispherical spectral emissive power moves 
to shorter wavelengths. The spectral emittance does not 
have a maximum at these shorter wavelengths, thus al- 
though Fig. 9 does not show this effect, the overall hemi- 
spherical emittance will eventually begin to decrease as 
the temperature increases. In the case of £ dependence on 
sheet thickness, T, e approaches an asymptote as t -» 


Conclusion 

The emittance of the thin liquid sheet in Fig. 1 was 
determined by two methods. The first method was to 
calculate it from transmittance data. The second method 
used the measured temperature drop across the sheet to 
determine emittance. For purposes of computing the emit- 
tance from the temperature drop, Tj - T 2 , the sheet was 
treated as a constant temperature device. It was shown 
that no significant error in the emittance results from this 
approximation. 

From the analysis and results, several points can be 
made. It is evident that the liquid sheet functions well as 
a radiator. The experimental emittance of the sheet is 
between 0.74 and 0.85, depending on the sheet thickness 
and oil temperature. A large emittance is needed for a 
space radiator, making the liquid sheet an ideal choice. 
There was also close agreement in the emittance value 
between the theoretical and experimental method. A con- 
cern in the beginning of the experiment was that measur- 
ing the small temperature drop would prove too difficult. 
Getting a very accurate measurement was necessary and 
temperature fluctuations in the oil could easily cause an 
incorrect measurement. By using thermistors with accu- 
racies of 0.01 K, the measurement was sufficiently accu- 
rate. To eliminate temperature (and velocity) fluctuations 
the oil flows into a large plenum before passing through 
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the slit. The plenum sufficiently dampens any fluctua- 
tions. An analysis of the fluctuations of the oil showed 
the fluctuations to be negligible or well over an order of 
magnitude less than the temperature drop due to radiation. 

The large emittances (e > 0.8) indicate that a liquid 
sheet makes an excellent low mass space radiator. Be- 
fore an actual liquid sheet radiator can be fabricated, 
however a collection method for the flow must be de- 
signed. Since the flow coalesces to a point this should 
simplify the collection process. 
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Figure 3. — Schematic of oil and window configuration 
for transmittance measurement. 
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Figure 5.— Transmittance of potassium bromide window. Figure 6.— Extinction coefficient for Dow Coming 705 

silicone oil. 
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Figure 7. — Schematic of test facility. Liquid sheet radiator 
experiment. 




Figure 9. — Comparison of theoretical and experimental 
total hemispherical emittance for Dow 705 oil. 
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Appendix A — Extinction Coefficient Derivation 


i 


* 


4 


At each interface in Fig. A-l, equations can be writ- 
ten for the energy flux (the incident radiation, qj, and the 
outgoing radiation, q 0 , for each side of the interface.) 
The primed variables refer to the radiation on the right 
side of the interfaces of Fig. A-l and the unprimed 
variables refer to the radiation on the left side of the 
interfaces of Fig. A-l. The general form of the four 
equations for interfaces b = 1, 2, 3 and 4 is 

< lob =r b ( lib + ( 1 - r b) C lib ( A_1 ) 

( iob = r b c iib + ( 1 - r b) c 4 ( A - 2 ) 

q ib = Vi^ob-i ( exce P l <Ia = *) ( A - 3) 


T = 




(i - rjr 2 T^ )(l - r 3 r 4*3 ) - x|(r 2 + fjtj [l - 2r 2 


3 + r 4 T 3 


[-*,]) 


(A -8) 


Because r^ = r 4 , r 2 = r 3 and Xj = X3, Eq. (A-8) can be 
simplified. If ri = r 4 = r aw (reflectance at air-window 
interface), r 2 = r3 = r mw (reflectance at media- window 
interface), and Xi = T3 = x w (internal transmission of 
window), then 



<4 = x b%b+i ( exce pt q| 4 = °) ( A - 4 ) 


(A -9) 


where q is the energy flux, or the energy per unit area per 
unit time, Tjis the transmission of window I, T2 is the 
transmission of the media, T3 is the transmission of win- 
dow II, ri and r 4 are the reflectances at the air-window 
interfaces, and r 2 and r3 are the reflectances at the media- 
window interfaces. 

These four flux equations are used to write energy 
balances for each of the four interfaces. The sixteen 
resulting equations can be written as a system of equa- 
tions. q^ 4 = T where T is the total transmission through 
the media to be tested and through the plates (windows) 
in which the media is contained. Algebraic manipula- 
tions show that, 

904=7^4 =rV T (A_5) 

1 ” r 4 1 “ r 4 


where x 2 has been replaced by x (transmission of media) 
in Eq. (A-8). Equation (A-9) is a quadratic equation for t 
as a function of T, r mw , r aw and x w , 

^(Vv +r aw X l[ 1 - 2r mw]) 2 +xx w( 1 - r aw) 2 ( 1 - r mw) 2 
-T(l-r a w x mw x2 w) 2 = 0 (A- 10) 

Thus, 

a = T(r mw +r aw x2[l-2r mw ]) 2 (A-ll) 

k = T w(l - r aw) 2 (l “ r m W ) 2 ( A - 12 > 


The system of equations is then solved forq^ and yields, 


zO-iX'-Ofr-O 


(' - Wl X' ‘ + vft - 2r 2 ])( r 3 + Vj[' " ^s]) 


Since, 



(A -6) 
(A -7) 


the transmission is solved for using Eqs. (A-6) and 
(A-7). This yields, 


: — T (! — r aw r mw X w ) 2 


(A -13) 


x = - 


— b±v/b 2 — 4ac 


2 a 


(A -14) 


Since, Tx = e _ad , then the extinction coefficient can be 
written 


. | — b ± Vb 2 -4ac 

I 2a 


a x=- 


(A -15) 
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(A -18) 


To solve for the extinction coefficient, the transmission 
of the window, x w , must also be calculated. It can be 
calculated using Eq. (A-8) illustrated above except in 
this case there are only two interfaces and the window is 
the new “media”. 

It can be seen that, T\ = = 0, r 2 = r3 = reflectance at air- 

media interface, T\ = T3 = 1 and t = transmission of 
window. The equation for the transmission of the win- 
dow becomes, 


T X ( } — r 2 )(l ~ r 3 ) 


1-T 2 r 2 r 3 


(A -16) 


Therefore, since r 2 = T 3 = r aw , 


a — T r 2 

a “ W aW * 

b = (l-r a w) 2 (A -19) 

and 

c = -T w . (A -20) 

The internal transmission of the window is written then 
as 




~ r awf 1 ~ r aw ) 4 + 4T v 


2 r 2 

w aw 


2T x L 

w aw 


(A -21) 


^Vaw+^-r a w) 2 -T w =0 (A -17) 

and 


The internal transmission of the window as a function of 
wavelength is used in Eqs. (A-14) and (A-15) and the 
extinction coefficient of the oil is then calculated and 
used to find the spectral emittance. 


Window I Media Window II 



^14 


© © © © 
Figure A-1 —Energy flux of absorbing media between 
transmitting windows. 
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Appendix B — Limit Calculations for Sheet Temperature Drop 


* 


V 


4 


In the limit of low background temperature (T M — > 0), 
the energy equation no longer has the solution (7). In- 
stead, we must solve the modified energy equation: 

^J = -2hTf 3 dz (B — 1) 

which has the explicit solution: 

^- = -(6hz + l)- 1/3 (B-2) 

A 1 

The solution Eq. (B-2) can be substituted directly into 
the power balance Eq. (1 1). This yields: 

/•6h+l 

xj [(6h + l)u -4/3 - u -1/3 ]du (B-3) 

The integration of Eq. (B-3) is straightforward and re- 
sults in 


AT 4h + l-(6h + l) 

im = - 


2/3 


lim 

T ^o hT l 


4h 


(B — 4) 


While the low sink temperature limit is appropriate 
for space applications, typical experimental conditions 
may have sink temperatures relatively close to the fluid 
temperature (such as the room temperature, for example). 
It has been found numerically that the percent error due 
to the constant temperature approximation has a mini- 
mum at T^ = 0 and increases slightly with increasing T M . 
Thus it is desirable to verify that this percent error never 
exceeds the limit of 8h/3 regardless of the sink tempera- 
ture. This can be done by considering the limit T*, — > Tp 
Two new parameters, <j> and %, are defined in the 
following way: 

^ = l-<t>, for $ -» 0,x -» 0 (B-8) 

A 1 A 1 

In the constant temperature approximation, <|> = 0 and 
Eq. (B-8) can be substituted into Eq. (1 1) to obtain: 

(B-9) 

Being interested only in limits, % is taken to be arbitrarily 
small and terms of higher than first order are neglected. 
The “exact” solution is found from the energy Eq. (12), 
which becomes: 


Since the range of h is 0.0004 to 0.004, or h«l, the 
temperature drop is found to be: 


dT/T 


d<j> 


(l-<f) 4 -(l-z ) 4 


= -2h dz 


(B-10) 


AT 4h + l-(6h + l) 2/3 , ^ ^ 

— = T5 — -1 (B-5) 


lim 
T -*> hT l 


4h 


The relative error in AT resulting from the constant tem- 
perature assumption (that is, AT = hTj) was found nu- 
merically to be proportional to h. In the limit: 


lim (AT>approx 


"(AT), 


exact 


(AT) 


approx 


_ 4h 2 -4h-l + (6h + l)^ _8 [i 
4h 2 3 

As h increases, this ratio decreases asymptotically to 
zero, so the following inequality holds: 


lim (AT)approx 


(AT) 


(AT)exact ^8h (B — 7) 


approx 


This equation has a closed-form solution, namely: 

<t> = x(l-e- 8h5 ) (B — 11) 

Substituting Eqs. (B-l 1) into (11) and eliminating higher 
order terms yields: 


lim *1 
T T, hT, 


l Jo 


— = 2 [(l-«l>) 4 -(l- X ) 4 ][l-z]dz^2 


f 4[-<!? + %][1 - z]dz = 8% f e 8hz [1 - z]dz 

JO JO 


8h 2 [8h-l + e _8h 


(B — 12) 
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Combining Eqs. (B-9) and (B-12), the limit on the error 
in making the constant temperature approximation is 
determined by: 


lim 

T T 


approx (AT) exact 

approx 



(B-14) 


lim 

To^T, 


^ AT ^approx ( AT ^exact 


approx 


32h 2 -8h + l-e 8h 8 

2 = “ h 

32h Z 3 

(B-13) 


Therefore, the results of Eqs. (B-7) and (B-14) hold for 
any < T\. 


As h increases, this ratio also decreases asymptotically 
to zero, so there is an upper bound equivalent to 
Eq. (B-7); that is, 
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